function [] = importCommutingMatrix2Matlab(dataFile, geoLabelsFile, outputFolder)
	
	% Add paths
	addpath('../code/general_funcs');
	
	% Read file with geographic labels
	opts = detectImportOptions(geoLabelsFile);
	opts = setvartype(opts, opts.VariableNames{1}, 'char');
	geoLabels = readtable(geoLabelsFile, opts);
	geoLabels = table2cell(geoLabels(:,1));
	
	K = length(geoLabels);
	
	% Read data file
	[mydata, ~] = makeArrayFromCSV(dataFile, 2, {'char', 'char', 'double'}, 0, {geoLabels, geoLabels});

	% Save un-normalized matrix
	W = mydata;   % NumGeoUnits (residence) x NumGeoUnits (travail)
	outputFile0 = sprintf('%s/MobilitePro.mat', outputFolder);
	save(outputFile0, 'W', 'geoLabels', '-v7.3');
	
	% Set diagonal equal to zero
	mydata(get_diagonal_idxes(K)) = 0;
	
	% Row normalize (and transpose)
	W1 = rowNormalize(mydata); % NumGeoUnits (residence) x NumGeoUnits (work)
	W2 = rowNormalize(mydata'); % NumGeoUnits (work) x NumGeoUnits (residence)
	
	% Save matrices
	outputFile1 = sprintf('%s/MobilitePro_Within_ResidenceLocation.mat', outputFolder);
	outputFile2 = sprintf('%s/MobilitePro_Within_TravailLocation.mat', outputFolder);
	
	W = W1;
	save(outputFile1, 'W', 'geoLabels', '-v7.3');
	
	W = W2;
	save(outputFile2, 'W', 'geoLabels', '-v7.3');
end
